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Distribution of work in isothermal non-equilibrium processes 

Thomas Speck and Udo Scifcrt 
//. Institut fiir Theoretische Physik, Universitdt Stuttgart, 70550 Stuttgart, Germany 

Diffusive motion in an externally driven potential is considered. It is shown that the distribution 
of work required to drive the system from an initial equilibrium state to another is Gaussian for 
slow but finite driving. Our result is obtained by projection method techniques exploiting a small 
parameter defined as the switching rate between the two states of the system. The exact solution 
for a simple model system shows that such an expansion may fail in higher orders, since the mean 
and the variance following from the exact distribution show non-analytic behavior. 
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In macroscopic thermodynamics, the work W spent in 
changing the state of a system at constant temperature 
T obeys 



W^ AF, 



(1) 



which is one version of the second law where AF is the 
difference in free energy of the final and the initial equi- 
librium state. As the system gets smaller, thermal fluc- 
tuations play an increasingly relevant role. Hence, this 
work acquires a stochastic contribution, i.e., the work 
follows a distribution function P{W). The shape of this 
function depends on how the system is driven. If this 
change is induced by the time variation of an external 
control parameter A(<), the distribution P{W) becomes 
a functional of A(<). 

Such distributions have recently become accessible ex- 
perimentally for systems with only a few degrees of free- 
dom diffusing in a thermal environment under the influ- 
ence of an externally controlled potential. Paradigmatic 
examples include dragging a coUoidparticle by an opti- 
cal tweezer through a viscous fluid jlj, I3 and the forced 
unfolding of RNA hairpins 13|. In both cases some re- 
alizations of the process show W < AF. Slightly over- 
stated, such findings have been called violations of the 
second law [l|. In a more conservative interpretation of 
the second law for such mesoscopic systems, the average 
work should and does still obey 



W = 



+ 00 



dW P{W)W ^ AF. 



(2) 



Obviously, the distribution P(W) is of paramount 
importance for a better understanding of isothermal 
stochastic dynamics. Exact statements, however, about 
P{W) are scarce. In 1997, Jarzynksi has shown un- 
der rather mild assumptions that the distribution P{W) 
obeys an integral constraint 



+ OC 



dW P{W)e 



-pw 



-PAF 



(3) 



for any external protocol A(t) |J,MjIS|- Here, /3 = l/fceT 
with Boltzmann's constant fee- Since this remarkable 



relation allows one to extract equilibrium free energy dif- 
ferences from measuring or calculating the work distribu- 
tion in non-equilibrium experiments or simulations, it has 
found widespread applications recently pi ItI IM liA llOl llll | . 
The statistical and convergence properties of this non- 
linear average deserve particular attention |l2l ll.'j . |lj| . 

For time-dependent quadratic potentials, i.e. linear 
stochastic equations of motion, P{W) can easily derived 
to be Gaussian [l5j. Jarzynski's relation then implies 
that the mean W and the variance a^ are necessarily 
related by [3| 



W = AF + l3(T^/2. 



(4) 



For these potentials, the Gaussian nature holds indepen- 
dently on the speed of the driving, i.e. independently of 
how far the system is from equilibrium. 

The purpose of this paper is to add a third general 
statement about P{W) to this list of exact results. We 
will show that this distribution becomes a Gaussian for 
slow but finite driving even if the equations of motion 
are non-linear. Since our approach is constructive, it 
yields an explicit algorithm of how to obtain the mean W 
and the variance a^ of this Gaussian distribution. In the 
quasi-static limit of infinitely slow external manipulation, 
this Gaussian reduces to P{W) = S{W - AF). 

A Gaussian character of the distribution P{W) near 
equilibrium seems to be expected or taken for granted in 
the recent literature 0, 0, llfl ll^ ll^ ■ Closer scrutiny 
of the references usually cited for this assumption - if 
any are cited at all -, however, reveals that they do not 
provide an explicit proof of this statement. The often 
cited papers by Hermans [13 and Wood et al. [13 ex- 
plicitly assume a Gaussian shape. Alluding in a more 
general way to an Onsager-Machlup functional [l9j also 
fails, since this Gaussian functional is derived for linear 
stochastic equations of motion. The presumably most 
promising case to date in favor of a Gaussian distribu- 
tion suggests to invoke the central limit theorem for the 
increments of work [20j. However, without translating 
this proposal into a definite calculation, which seems to 
be non-trivial for time-dependent potentials, this argu- 
ment is not a clear-cut proof yet, let alone does it give an 
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expression for mean and variance of this putative Gaus- 
sian. 

Based on this unsatisfying state of affairs concerning 
such a fundamental issue, we beheve that a constructive 
derivation of the Gaussian nature of this distribution for 
finite but slow driving is indeed called for as a step to- 
wards a comprehensive theory of isothermal stochastic 
dynamics. 

For the derivation we consider a finite classical sys- 
tem coupled to a heat reservoir of constant temperature. 
Let then x = (xq, . . . , a;„) be the state of the system with 
energy V\ (x) where A is an externally controlled parame- 
ter. The stochastic dynamics is governed by the Langevin 
equations |21| 



-My 



dx^ 



n.{t) 



(5) 



where /i^ are the mobility coefficients |2J| and rii{t) is 
the thermal noise representing the heat bath with 

{r,,{t)) = and (ry,(t)7?,(i')) = '^^i^At - t'), (6) 

where (. . .) denotes the ensemble average. We describe 
the continuous process of switching the system from an 
initial state (A(i = 0) = 0) to a final state (A(i = ig) = 1) 
by a protocol X{t), over a total switching time ig- With- 
out loss of generality, we set 



A(t) = t/U 



(7) 



and hence a constant switching rate X = t^^ [j|. 

We now consider an ensemble of infinitely many real- 
izations of this Markov process, each evolving stochasti- 
cally according to eq. Q- The normalized distribution 
of this ensemble in phase space /(x, t) obeys a Fokker- 
Planck equation [21| 



dtf ^ Lxf with Lx = ^My 



dVx 
dx-i 



1 d 
P dxj 



(8) 



equivalent to the Langevin equations ^ . This introduces 
the (through A time-dependent) Fokker-Planck operator 
Lx- The stationary solution of eq. (jSJ for fixed A is the 
equilibrium distribution 



/r(x) 



-PVxi^)/ / j^/ g-/3yA(x')^ 



(9) 



The total work performed along one particular trajec- 
tory x(i) up to time t is the time integral 0, |22| 

W^[x(t),i]^£di'A^(x(i')). (10) 

We can now compose a combined stochastic process con- 
sisting of {x, W} as [13 

dVx 



W =A 



dVx 



dX 



(11) 
(12) 



Note that the equation of motion for W does not have 
an independent noise but is stochastic through the x- 
dependence of Vx- The joint probability distribution 
function p(x, W, t) then obeys a Fokker-Planck equation 



dtp^ Lx + xlY p 



where 



fw 
^x 



dVx d 
dX dW 



(13) 



(14) 



represents a drift term of the work. Eq. (|13|l is exact, no 
matter how far the system is driven out of equilibrium. 
The reduced probability distribution of the work P{W, t) 
can be obtained by integrating out x as 



P{W,t) = f dxp{x,W,t). 



(15) 



Since we start the process out of thermal equilibrium, 
the x are initially distributed according to the canonical 
distribution and therefore the initial condition is 



pi^,w,o) = r,\^)s{w) 



(16) 



As our main theoretical tool we introduce a projector 
11;^ acting on a function 0(x, W, t) such that 



nA0: 



fcq 
JX 



dxXx',VK,t). 



Note that 



LxTix<t> = ^xLx 



0. 



(17) 



(18) 



The first statement is evident from definition H17|l and 
the fact that f'^ is in the null space of Lx- The sec- 
ond conclusion follows when <f> is expanded in terms of 
eigenfunctions of Lx- Then the Fokker-Planck operator 
Lx cancels the eigenfunction to eigenvalue 0, which in 
fact is f'^, whereas the projector annihilates all other 
eigenfunctions corresponding to higher eigenvalues. 

The other important property of the projector Ha, 
which distinguishes it from the usual app lication to the 
adiabatic elimination of fast variables [2l| , is that it does 
not commute with the time derivative but rather leads 
to 

[duiix] = dtilx ' ii-xdt = -XpSxii-x, (19) 

where we define 

dVx IdVx 



Sx — 



dX \ dX 



(20) 



The equilibrium ensemble average {■ - -) x i^ defined as 

(0),^y'dx/^<^(x)<^(x,i). (21) 



We can now expand the joint probability p(x, W, t) 
for small A, which corresponds to a separation of time 
scales |23|. The slow time scale is A = i/ig. The fast time 
scale, which we do not need explicitly, is determined by 
the intrinsic relaxation processes. As the time derivative 
transforms according to dt -^ Xd\, switching to the slow 
time scale eq. (|13|l becomes 



dxP 



x-'l^ + lY p. 



(22) 



By using the projector !!> we decompose the distribution 
function p = po + pi into 



and 



Po(x, W, A) = lixP = fr{^)P{W, A) (23) 



Pi{^,W,X) = {l-Ilx)p. (24) 



We apply 11^, respectively (1 — Ha), to eq. H22|l and keep 
in mind both cq. (|18|l and the commutator H19|) . We 
finally get the two coupled differential equations 



dxpa = ^aPo + ^aPi - /^-S'aPo 



(25) 



dxpi 



Pi 



X-'Lx + Ai 

where we abbreviate 

riAir and Al = {l-tlx)L 



AIpo + (3Sxpo (26) 



Al 



(27) 



In this form, an expansion in A becomes possible. In 
lowest order (A — > 0), eq. (^5)) imphes XaPi = 0. Since pi 
is orthogonal to the null space of Lx by definition (|24|l . 
Pi = follows. For a solution of cq. (P5|l wc explicitly 
calculate 



dVy 



d 



^A/r = -/r dx ^/r^^ = -ir i^) i^ 



dX ■"- dW 



Using this and eq. (|23|l we finally obtain 



dVx 
dX 



d 



dP 



dVx 
dX 



dP 
dW 



dW 

(28) 



(29) 



for the distribution of the work P{W,X). The solution 
of this equation is P{W) = S{W — AF) with the initial 
condition P{W, 0) — 6{W) following from eq. H16|l . where 



we recognize 



AF 



. - m\ 



(30) 



as the change in free energy of the entire process. We 
thus have recovered for A ^ the quasi-static limit as 
expected. 

To first order in A we get from eq. H26|l 



Pi = -ALi Al + f3Sx Po 



(31) 



Putting this back into eq. H25|l we get after a straightfor- 
ward calculation a diffusion type equation for P{W, A) in 
the form 



dP 
9A 

where 



dVx 
dX 



XpSx 



dP 

dW 



XS- 



d^P 



Sx 



dx ^L-^Sxfx-. 



(32) 



(33) 



The solution is a Gaussian 



P{W) 


1 


2a2 




(34) 


V2w2 P 


with variance 








a^ = 2X f dX Sx 
Jo 


(35) 


and mean 






W= / dA 
Jo 


mh^^^^] 


= AF+|a^ 


(36) 



This is the central result of the present paper [2^. First, 
it proves that the distribution of the work in isother- 
mal non-equilibrium processes is Gaussian in the near- 
equilibrium regime. Second, we recover independently 
from Jarzynski's relation (PJ the constraint that the mean 
and variance are connected according to eq. Q). Third, 
it yields an explicit algorithm of how to calculate these 
quantities. 

For an assessment of the range of validity of this ap- 
proximation, we recall that it is based essentially on a 
separation of time scales. Hence, the Gaussian distribu- 
tion will be a good approximation as long as 



At < 1, 



(37) 



where r is an intrinsic relaxation time. 

As an example, we illustrate our approach for a simple 
one-dimensional case, where we can compare our expan- 
sion with an exact solution [l5j . We consider a colloidal 
particle at position x with mobility ^ trapped by an opti- 
cal tweezer whose center y{X) is moved at constant speed 
V through a viscous fluid (see FIG.^. The potential of 
the trap is assumed to be harmonic near the focal point 



Vx{x) = {k/2){x-y{X)f 



(38) 



with effective strength k. In this case, the free energy is 
independent of y{X). For the two states, we choose with 
y{X = 0) = and y{X = I) = L two positions of the 
trap. The switching rate becomes A ~ v/L, while the 
relaxation time is r = 1/^fc. 

Within our scheme, wc first have to calculate ^a = 
I/A Sxfx^ in cq. H33() which amounts to solving the inho- 
mogeneous differential equation 



La^> 



sxir. 



(39) 




FIG. 1; Scheme of the experiment of Ref. [J]. A colloidal 
particle at position x{t) is dragged by an optical tweezer with 
focus at y{X) = LX through a viscous fluid. The effective 
potential is harmonic with a spring constant k. 



With the Fokkcr-Planck operator from cq. © we get ex- 
phcitly 



d_ 

dx 



'^(^-^^) + ^^. 



*A = 



L(x — La) exp 

27r 



--I3k{x- LXf 



(40) 



where A only appears as a parameter. This is easily 
solved as "^ \{x) = f'^{x)Lx/ j-i and thus the average work 
becomes 



W = L^k At. 



(41) 



Of course, for this harmonic potential, the distribution 
P{W) is Gaussian at any driving [l^. The exact result 
for the mean W^ as a function of A reads 



W = L^k 



Xt - AV^ ( 1 



-l/Ar 



(42) 



which agrees to first order in A with eq. (|41|) as expected. 

The exact expression (|42|l points to an interesting 
property which seems not to have been discussed yet 
in the context of stochastic dynamics |2g. The expo- 
nentially small last term shows that the average work is 
non-analytic in A. We expect that if our expansion of 
eq. H26() was extended to the next order, some signature 
of this non-analyticity should show up. Therefore the ap- 
proach to equilibrium (or the deviation from equilibrium) 
even in this almost trivial case is somewhat subtle. 

In summary, we have shown for general diffusive sys- 
tems that the distribution of work required to drive the 
system from an initial equilibrium state to another is a 
Gaussian for slow but finite driving. Its mean and vari- 
ance can be obtained from solving an inhomogeneous dif- 
ferential equation involving the Fokker-Planck operator. 
As an exactly solvable case shows, these quantities are 



non-analytic in the switching rate. This result indicates 
that in general calculating the next order correction to 
the Gaussian derived here may face fundamental difficul- 
ties. 

Stimulating discussions with O. Braun and R. Finken 
are gratefully acknowledged as well as valuable hints and 
comments by H. Spohn. 
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